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We study the relative motion of nearby free test particles in cosmological spacetimes, such as 
the FLRW and LTB models. In particular, the influence of spatial inhomogeneities on local tidal 
accelerations is investigated. The implications of our results for the dynamics of the solar system 
are briefly discussed. That is, on the basis of the models studied in this paper, we estimate the tidal 
I influence of the cosmic gravitational field on the orbit of the Earth around the Sun and show that 

, the corresponding temporal rate of variation of the astronomical unit is negligibly small. 

1—5 

a^ 

I. INTRODUCTION 

CJ'. 

The gravitational influence of distant galaxies extends over all of the bodies in our cosmic neighbourhood; therefore, 
5^ ' the relative motion of these bodies would be affected by the tidal field of distant masses. In general relativity, it is 
bJ), natural to study the relative motion of a test particle with respect to an observer in the quasi- inert ial Fermi co-ordinate 
system [1] established about the worldline of the observer. To simplify matters, let us assume that the observer follows 
' a geodesic worldline C and let ^'^{a) be the observer's local orthonormal tetrad frame that is parallel transported along 
^ ■ C. Here A''(o) = dx^'/dr is the observer's local temporal axis and r is its proper time, while i = 1,2,3, are 

(Nl ' unit gyro directions that constitute the observer's local spatial frame. We choose units such that c — 1 throughout 
T-H , this paper. An event P in the neighbourhood of the observer has Fermi co-ordinates X'^ = {T,X); in fact, P can 
be orthogonally connected to C at Pq via a unique spacelike geodesic of proper length a such that cr = at Pq. Let 
the proper time along C at Pq be r and = (da;^/dcr)Q be the unit tangent vector to the spacelike geodesic at Pq; 
' then the Fermi co-ordinates of the event P are defined by T = t and — a A^(*^ It follows that the observer 
, permanently occupies the spatial origin of this Fermi co-ordinate system. 
The spacetime metric in Fermi co-ordinates can be expressed as 



""500 = -1 - ^Ro^OJ{T) X'X^ + ■■■ , (1) 
^9^^ = -^^Po,,fe(T)XJX'= + --. , (2) 

^ ■ ^9^J = '5,, - Y^^k.i {T)X^X' + --- , (3) 

where 

^RalS-ysiT) = P^iypcr A^(q) A'^f^) A''(^) A'^(5) (4) 

is the Riemann curvature tensor along C projected on the tetrad frame of the observer. The Fermi co-ordinates are 
admissible in a cylindrical spacetime region around C with \X\ < C, where £ is a measure of the radius of curvature 
of spacetime. 

The motion of a free test particle is given by 

d^X^ p „ dX" dX^^ 



Kg-r--7- = 0, (5) 



ds2 ~^ "/^ ds ds 
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where dX^/ds = F (1, V) is the particle's timeUke four- velocity vector. Equation (5) can be expressed as 



1 dr dX" dXf' 

and the reduced geodesic equation 



r dT "'^ dT dT ^ ' 



This latter equation, to linear order in distance away from C, is the generalized Jacobi equation [2] 

+^ (3 ^Rokjo V + ^Rikji + ^Rokji V V^) X^ = 0, (8) 
and the modified Lorentz factor is given by 

= l-V^ + ^R,nojX'X^ + ^^Roj,kX^V' X^ 

+ \^RikjiV'X^V^ XK (9) 

Equation (8) reduces to the Jacobi equation when the velocity-dependent terms are nc;gUgibly small. 

In some situations of physical interest, it is possible to have purely one-dimensional motion, say along the Z- 
direction. In this case, equation (8) reduces to 

^+k{T){l-2Z^)Z = 0, (10) 

where Z = dZ/dT and k{T) = ^ Rtztz- For \Z\ <C 1, equation (10) reduces to the standard Jacobi equation. 
However, when the relative speed cannot be neglected in comparison with the speed of light, the generalized Jacobi 
equation (10) has solutions with Z = ±l/\/2 such that the relative motion within the linear approximation scheme 
is uniform at the critical speed Vc = Below this speed, the relative motion should essentially conform to 

expectations based on the post-Newtonian approximation. On the other hand, novel relativistic tidal effects are 
expected to occur above the critical speed. 

Let the critical solutions of equation (10) be expressed as Zc{T) = Zi ± Vc{T — Ti), where Ti is an initial time and 
Zi = Zc{Ti). Then, the behaviour of Z{T) near a critical solution can be examined by introducing C(^^) = Z — Z^- It 
follows from (10) that 

|| - 2 k{T) (C^ ± 2V, C) (C + ^c) = 0. (11) 
Keeping only terms of linear order in (11), we find 

\% = ±4K^cfc(r), (12) 

where C, = d^/dT. It follows that the long-term behaviour of this motion can be determined from 

iT+Co)k{T)dT, (13) 

where Co = — Tj ± V^Zi and ±2VcZc = T + Co- For instance, the critical solutions are attractors for X = — oo [2]. 
Equation (10) has been extensively studied in black hole spacetimes in connection with the problem of astrophysical 
jets [3, 4]. 

The main purpose of this paper is to study equation (10) for some cosmological models. It turns out that for the 
standard FLRW models, k{T) is characterized by the deceleration of the universe. This point is demonstrated in the 
following section. In section III, a locally inhomogencous model is employed to determine the corresponding k{T). 
In this way, the influence of cosmic spatial inhomogeneities on local measurements is elucidated. This treatment is 
extended in section IV to the Lemaitre-Tolman-Bondi (LTB) models. The consequences of the results of section IV for 
the dynamics of the solar system are briefly discussed in section V and Appendix A. Section VI contains a discussion 
of our results. 
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II. FLRW MODELS 



Let us first consider the FLRW spacetimes given in isotropic co-ordinates = {t,x^) by 

P{r) 



ds^ = -dt^ + — f4 % da;Ma;^ (14) 



where S{t) > is the scale factor and 

/(r) = l + \i^r^ (15) 

with K = — 1, 0, or +1 for the open, flat, or closed universe models, respectively. We note that the spatial co-ordinates 

a;' are dimcnsionlcss, while the scale factor S has the dimension of length. 

A fundamental observer in this spacetime occupies a fixed position in space and follows a geodesic worldline C with 
proper time t = t. Let us choose a fundamental observer with A'*(o) = '^'^o and ^^(i) = [f{'r')/S{t)] 6^i, which can be 
shown to be parallel transported along the geodesic worldline. The curvature tensor ^ Rap-yS along C can be expressed 
as a 6 X 6 matrix {^TZjk), where the indices J and K range over the set {01, 02, 03, 23, 31, 12}. Then, 

where for the FLRW models E = u{t) I, B = 0, and A'' = w{t) I. Here / is the 3x3 identity matrix and 

, , 1 d^S „, 

^0 = -^^, (17) 

1 (dSV K 



We note in passing that (17) and (18) imply 

^ dw „ , , d5 

5-=-2(. + .)-. (19) 

It is possible to show that the only nonzero components of the Einstein tensor, G^i, = — ^ g^j, R, are given by 

G^.A^(o)A'^(o) = 3 wit), (20) 
G^.A^ = [2u{t)~w{t)]5,,. (21) 

The gravitational field equations (with a cosmological constant A) are given by 

G^.+Ag^, = SttGV, (22) 

where T^i, is the stress-energy tensor for a perfect fluid 

Tij,i, = Kt)ui^u^+p{t){gfj,^ + Ui^u^). (23) 

Here, /^(t) and p{t) are the invariant density and pressure of the cosmic fluid, and the fundamental observers are 
assumed to be comoving with the cosmic fluid, i.e. u'^ = X^{o)- Equations (20)-(23) imply that u and w, deflned 
respectively by (17) and (18), are given by 

AttG , „ , a . 
u{t) = — (M + 3p)--, (24) 

, , 87rG A . 
w{t) = — (25) 

It follows from (19), (24), (25) and the constancy of A that 

1 du. 1 dS , , 
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This is the expression of the first law of thermodynamics for the adiabatic flow of the perfect fluid under consideration 
here. The standard Hubble and deceleration parameters {H and q) are given by 

so that u = qH"^ and w = H"^ + k/S^. 

A Fermi coordinate system (T, X) can be established along the worldline of a fundamental observer; then, the 
observer has Fermi coordinates (T, 0), where T = t. The FLRW spacetime metric in spherical Fermi co-ordinates 

X = p sin G cos $ , F = p sin © sin $ , Z = p cos 6 , (28) 



can be expressed as [5] 

-ds^ = - [1 + u{T)p'^] dT^ + dp2 + 



p2(de^+sin2ed$2). (29) 



This is based on (l)-(3) and thus holds only to second order in distance away from the observer. It appears that 
explicit Fermi co-ordinates can be constructed for FLRW models based on the results of Ref. [4]; however, the scale 
factor S{t) must then be explicitly specified. To have a general treatment for arbitrary S{t), we must use the implicit 
approach as in (l)-(3). 

The spatial isotropy of the FLRW models implies that an equation of motion of the form of equation (10) is possible 
along any direction in space with fc(T) = u{T) = qH^. Thus the nature of such motion is characterized by the sign 
of the deceleration parameter. Imagine, for instance, a free test particle moving radially away from the observer. For 
an initially infrarelativistic particle with recession speed less than Vc, the particle decelerates (accelerates) for g > 
{q < 0), as expected. However, for an initially ultrarelativistic particle with recession speed above Vc, the particle 
accelerates (decelerates) for > (g' < 0). In the following sections, we simply concentrate on the nature of k{T) in 
certain inhomogeneous models. 

In the slow-motion approximation, it is possible to extract a "Newtonian" gravitational potential [5] 

Vn = ^qH^P^ (30) 

from (29). Equation (30) has been employed in the discussion of the influence of cosmology on local phenomena [6]. 
Ref. [6] contains only a partial list of papers on this subject; for background material and further references, see [7]. 
The form of the quadratic potential (30) suggests that the influence of the cosmic gravitational field is reflected in 
local experiments in the solar system, for instance — via a relative tidal acceleration of the form ^cosmos = —qH^X. 
It is interesting to note that the deceleration parameter q also appears in the second-order expansion of luminosity 
distance rf^ versus redshift z\ that is 

dL = ^z+^{l-q)z^ + --- . (31) 

This "degeneracy" is removed by spatial inhomogeneities as demonstrated in the following section. 

It appears highly likely that instead of providing evidence for dark energy, the observational data from type la 
supernovae have demonstrated the inadequacy of the standard spatially homogeneous FLRW models of the universe. 
That is, spatial inhomogeneities must be taken into account when comparing observational data with theoretical 
models of the universe. This is based on the fact that spatial inhomogeneities mimic dark energy [8]. A recent useful 
review of this topic is contained in [9] . 



III. INHOMOGENEOUS MODELS 



It is important to investigate how spatial inhomogeneities in the universe are reflected in local spacetime measure- 
ments. We begin with a metric of the form [5, 10] 

-ds^ = -a2(^,r)d^2 + 62(^^r)dr•2+7^2(^,r)(d02^sin2 0d</.2), (32) 
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where 



a{t,r) 
b{t,r) 

n{t,r) 



S{t), 
Sit) 



(33) 

(34) 

(35) 



Here, a, /3, and 7 are parameters in the expansion of the cosmic gravitational potentials in powers of the dimensionless 
radial co-ordinate r away from the centre of spherical symmetry (r = 0). This model has been discussed at length 
[5, 10]; in particular, it has been shown in [10] that this local inhomogcncous model is generally valid for a cosmic 
fluid that has pressure and satisfies a reasonable equation of state of the form p = p{ii) with p > 0, /x > and ii> 3p 
OTc n> p. We note that at the order of approximation indicated in (33)-(35), the homogeneous models of the previous 
section are recovered for a = and /3 = 7 = — k/2. 

We imagine, as before, observers that are fixed in space and carry orthonormal tetrad frames along their worldlines. 
The natural tetrad of such an observer is diagonal with nonzero elements A'^(o) = a~^, X^m = b~^, = T^~^, and 



AS 



(3) 



{TZ sin I 



Wc arc interested in a Fermi co-ordinate system along the worldlinc C of the observer that is at 



the centre of spherical symmetry. It can be shown explicitly using equations (32)-(35) that C is a geodesic and the 
spatial frame is indeed parallel propagated along C in the r — )• limit. Alternatively, one can employ local Cartesian 

co-ordinates as in the previous section. In any case, the curvature components ^RapfS along C can be explicitly 
calculated and the results can be expressed as a 6 x 6 matrix (16) with 



E = U{t)I, 
where for this inhomogeneous model 



B = 0, 



N = W{t)I, 



U{t) = qH 



W{t) = + 



2 , aft) 
^2 ' 

2 , Pit) - 37(i) 
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(36) 

(37) 
(38) 



Here, q and H are defined in terms of the scale factor Sit) as in (27). Moreover, the analogues of equations (24) and 
(25) hold for fiit, r) and pit, r) evaluated at r = 0; that is. 



Uit) = ^[Mi,0) + 3p(t,0)]-|, 

Wit) = ^MM) + |. 



(39) 
(40) 



The spacetime metric in spherical Fermi co-ordinates about C has the same form as equation (29), except that 
uit) Uit) and w{t) — ^ Wit); moreover, t = T, since r = along C. Let us note in this connection that the quantity 
/3 — 37 in (38) is related to the spatial curvature and reduces to k in the homogeneous limit [5, 10]. The motion 
of a free test particle in the neighbourhood of C can be studied as in the previous section except that in this case 
A;(r) = J7(T) in equation (10). Furthermore, the influence of the cosmic gravitational fleld on local physics in the 
slow-motion approximation appears as a "Newtonian" acceleration of the form 



-H^ iq-C)X, 



where C it) is the inhomogeneity parameter given by 



C = -- 



iSH)' 



(41) 



(42) 



It is interesting to note that in this inhomogeneous model, the luminosity distance-redshift relation can be expressed 
as [5, 10] 



dr 



j^Z+^[l-iq + C)]z' + 



(43) 
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so that the measured "deceleration" parameter is in fact q + C. It follows that both q and C could be determined 
from observation if it were possible to measure k{T) = H'^{q — C) as well. It follows from (39) that 



H'^{q-C) = A-kG 



■M(i,0)+p(t,0) 



A 
3" 



(44) 



Hence, in the absence of a cosmological constant we have C < q, which is a noteworthy inequality. 
Finally, let us note that by imposing the comoving co-ordinate condition [10] 



dt 



(/3-37) 



-2H{t) a{t) , 



one can show from (39) and (40) that 



(45) 



d_ 
dt 



(46) 



which is the analogue of (26) in this inhomogeneous case. 



IV. LTB MODELS 

Imagine a spacetime metric of the form (32) with 

a{t,r) = 1, b{t,r) = / ^ ' ' , (47) 

where TV = dlZ/dr. It turns out that the gravitational field equations (22) are satisfied in this case for pure dust 

T^u = n{t,r)u^Uu (48) 
in comoving co-ordinates (i.e. = 6^0) provided 

Sir) = \v}-^--,n\ (49) 
^ ^ 2 7^ 6 ^ ' 

= 47^/.(^,r)7^2 7^^ (50) 

dr 

Here, £{r) with E > —1/2 has the interpretation of the net energy per unit mass of the spherical shell of dust at radius 

r, M{r) has the interpretation of mass within a sphere of radius r, and TZ = dTZ/dt. These spherically symmetric 
inhomogeneous dust models were first discovered by Lemaitre [11] and further studied by Tolman [12] and Bondi [13]; 
see [14] for a detailed discussion. 

It is a general result that in a spacetime with a metric of the form —dt^ + gij{t, x'') dx"^ dx^ , any test particle 
that is at rest in space follows a geodesic [2]. Therefore, the fundamental comoving observers in the LTB model 
follow geodesies. Moreover, each fundamental observer carries the standard orthonormal tetrad frame (discussed in 
the previous section) that has its axes along the directions of the co-ordinates employed in (32). It can be shown 
explicitly that such a tetrad frame is parallel transported along the geodesic path of a fundamental observer. The 
off-center cosmological measurements of such observers have been discussed in [15]. We are interested, however, in 
the local tidal dynamics of nearby test particles. 

It proves convenient, for the sake of simplicity, to orient the spatial axes of the Fermi frame such that the X-axis 
points along the polar (9) direction, the Y-axis points along the azimuthal (0) direction, and the Z-axis points along 
the radial (r) direction; that is, for the purposes of this section we choose A^(o) = = 5^o, A^(i) = (0, 0, 1/7?., 0), 
A''(2) = (0,0,0, 1/(7?. sin^)) and A'' (3) = (0, 1/6,0,0). We find that for a fundamental observer at a fixed position in 
space, the nonzero components of the curvature tensor can be obtained from 

^-Roioi = ■'^-Ro202 = Ki{T) , (51) 
^7?0303 = K2{T), (52) 

^i?3131 = ^i?3232 = i^3(r), (53) 

^i?i2i2 = Ki{T). (54) 
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Thus the analogue of (16) in this case is given by E = dia.g{Ki,Ki,K2), B = 0, and N = dia.g{K3, K3, K4). Here 

2 



n 
n 



Ka{T) 



K2{T) = Ks{T) 



nn' dr \2 



(55) 
(56) 



where t = T and r is simply a constant in the final expressions. Then (l)-(3) together with (28) imply that the Fermi 
system (T, X) along the worldline of any fundamental observer has a metric of the form 



[1 + (Ki sin^ Q + K2 cos^ e) p2] dr^ + dp2 



l--K3p']p'dQ' 



+ 



1--{K4^ sin^ e + Ks cos^ 6) 



sin2ed$^ 



(57) 



which is axially symmetric about the Z-axis (i.e. the radial r-direction) , as expected. Local dynamics in this spacetime 
is considered in Appendix A. 

We are particularly interested in the form of this metric in the neighbourhood of the fundamental observer at the 
center of spherical symmetry (r = 0). Assuming that the spacetime manifold is smooth at r = [13], we can write 



n{t,r) = rS{t) 



l + \^{t)T+\ll{t)r^ + 0{r^) 

I D 



Mir) = iM"'(0)r3 + O(r4). 
6 



Here, the scale factor is given by 5(t) = 7?.'(f, 0) > and 

- 7^'(^,o) ' 

Using these approximate expressions, we find that 



Kx = UiT) 

= >V(T) [1 + 0(r)] , 



1 + -T]r + 0{r^) 



m 

K2 
K4 



n"'it, 0) 
n'{t,o) ■ 



UiT) [l+r]r + 0{r 
W(r) [1 + 0(r)] , 



where 



S 

's 



w = - 



S"{0) 



(58) 

(59) 
(60) 

(61) 

(62) 
(63) 

(64) 



and rj = ^2S'A + SA^ / S. The Fermi metric for r = then takes the form 

-ds^ = -{l+Up"^) dT^ + dp^+(l-^yV p^^ p^ (de^ + sin^ 9 d^^) . 



(65) 



It is now a simple matter to recognize that the gravitational field equations for metric (65) result in the analogues of 
equations (39) and (40) in this pressure-free case, namely. 



47rG , A 
= — M(t,0)--, 



Here iJ.{t, 0) can be expressed as 



/x(t,0) = 



M'"{0) 1 
Stt S^{t) 



(66) 
(67) 

(68) 
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by virtue of equation (50). 

It follows from a comparison of (64) with (18) that —£"(0) plays the role of spatial curvature in these inhomogeneous 
models. Moreover, k{T) — qH^ in this case and it is interesting to compare this appearance of the deceleration 
parameter, defined in terms of S{t) in the standard manner, with the effective deceleration parameter Q obtained 
from the luminosity distance-redshift relation. It has been shown [10] that in the general pressure-free case, 



dL = 



where the inhomogeneity parameter is given by 



Q-q = 



1 dA{t) 



(69) 



(70) 



The LTB models have recently received much attention as simple alternatives to the dark-energy models of the 
universe; see, for example, [9, 16-20] and references therein. 
Equations (64) and (66)-(68) may be written as 



1 



1 /ds;y _ 1 

2\dt J ~ 6 
£"{0) 



GM"'{0) 



AS' 



H 



S' 



A 



(71) 
(72) 



The scale factor S{t) can be simply determined from (71) by quadratures. In general, however, a simple analytic 
solution is not available. The qualitative behaviour of S{t) can be studied using the interpretation of (71) in terms 

of the radial motion of a "particle" of unit mass with kinetic energy 5^/2 and total energy £"{0)/2. That is, a graph 
of the potential energy in (71) versus S can be used to illustrate the fact that motion can take place only in regions 
where 



£"{0) > - 



GM"'{0) 



+ AS' 



(73) 



Alternatively, (71) can be integrated numerically with the initial condition that at the present epoch to, S{to) = Rq, 
where i?o is the (spatial) curvature radius. Let us note that for £"{0) = and A = 0, 



S = 



-GM"'(0) 

4 ^ ' 



1/3 



(74) 



so that H = 2/{3t), and q = 1/2 follows immediately from (72); therefore, this case corresponds to the familiar 
Einstein-de Sitter model. Equation (72) is used in the next section to estimate the influence of the cosmic gravitational 
field on the dynamics of the solar system. 



LOCAL DYNAMICS 



The expanding universe is over ten billion years old; therefore, the cosmic tidal acceleration within the solar system 
is expected to be relatively very small. The long-term cosmological evolution of an "isolated" gravitationally bound 
system in an expanding universe is beyond the scope of this work; instead, we are interested in observable cosmological 
perturbations on the dynamics of the solar system. As explained in detail in the Appendix, the influence of spatial 
inhomogeneities on a Keplerian binary system can be studied on the basis of the Fermi coordinate system (57) 
associated with the LTB model. However, the small anisotropy of the cosmic microwave background radiation implies 
that a useful estimate of the effect can be obtained by ignoring any deviations from spherical symmetry. Therefore, 
we consider the quasi-inertial Fermi coordinate system that can be established along the worldline of the fundamental 
observer at r = 0. Within this co-ordinate system associated with equation (65), imagine the approximately elliptical 
relative orbit of a Keplerian binary system that is perturbed by the cosmic tidal acceleration 



d^X GMo 

dr2" ^ p3 



X 



(75) 



where at the present epoch, ^cosmos = ~<1oHq X according to equation (65). Here, Mq is the net inertial mass of 
the binary system. Let us note that if the members of the Keplerian binary were test particles, i.e. Mq 0, then 
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equation (75) would simply reduce to the Jacobi equation at the present epoch, since the relative speed of the particles 
is assumed to be negligible compared to the speed of light. If the external tidal acceleration is turned off at any instant 
of time;, the resulting orbit is the osculating ellipse at that instant with eccentricity e and semimajor axis A. Thus, 
the Newtonian orbital energy of the osculating ellipse is —GMq/{2A). On the other hand, it follows from (75) that 
the rate of change of this orbital energy is given by ^cosmos • so that 

GMq dA ^ 

"cosmos ^ ' V ' ^/ 

For the osculating ellipse, 

X-V = c.KA^evT^ ^/^'^ , (77) 

1 + e cos 

1 /2 

where ujk = (GMq/A^) is the Keplerian frequency, and tp is an azimuthal angle of the orbit ("true anomaly"). 
Moreover, it follows from (72) that 2qQ = 1 — where xo is the present value of 



1 

X 



(78) 



Putting equations (76)-(78) together, we find that AA/dT has the same periodicity as the orbit and its average over 
an orbit vanishes; that is 



d^ HqAb,^ , Vl-e^ siny 

dT = 1 + ecosy • (''^ 

This is in fact the Lagrange planetary equation for the semimajor axis of the osculating ellipse in this particular case. 
It turns out that — on the average — the cosmological perturbations under consideration here cause a precession of the 
orbit in its plane, but leave the orbit otherwise unchanged. The Appendix should be consulted for further details 
about the average behaviour of the orbit in this case. 
For the orbit of the Earth around the Sun, 



H^Ae 1 



Wk 4 



xlO-^cm/yr, (80) 



where e « 0.02, A « 1.5 x 10^"^ cm, and Hq « 70 km s~^ Mpc~^. Even if — as a result of spatial inhomogeneities 
or a nonzero cosmological constant — the absolute magnitude of the deceleration parameter is enhanced by several 
orders of magnitude or so, the rate of variation of the astronomical unit would still be too small to be detectable 
at present. This should be contrasted with the recent reported secular increase of the astronomical unit, based on 
radiometric data, amounting to about 10 cm per year [21, 22]. We therefore conclude that the presence of cosmological 
inhomogeneities (or a cosmological constant) cannot change the conclusion that the expansion of the universe has a 
negligible influence on the dynamics of the solar system. 

The results of this section as well as Appendix A indicate that solar-system anomalies [22] such as the Pioneer 
anomaly — cannot be explained in terms of cosmological perturbations based upon the general relativistic cosmological 
models considered in the present work. This is particularly evident from the negligibly small magnitude of ^cosmos 
over the solar system, since for |X| - 100 AU, Hi \X\ - 10"^° cm/s^. This is about thirteen orders of magnitude 
smaller than the anomalous acceleration of Pioneer spacecraft. 



VI. DISCUSSION 



For local systems, such as the solar system, each body is subject to the gravitational influence of the whole mass- 
energy content of the universe; therefore, the relative motion of bodies is only affected by the tidal acceleration of the 
cosmic gravitational field. The main purpose of this paper has been to study the general features of tidal dynamics in 
cosmological models. Particular emphasis has been placed on inhomogeneous models, since spatial inhomogeneities 
mimic dark energy. Indeed, we have elucidated the influence of inhomogeneities on tidal dynamics. The results of 
this work could therefore be of interest in the theoretical study of the tidal interaction between galaxies. 

The tidal influence of the cosmic gravitational field on the solar system can be estimated and the result turns out 
to be too small to be measurable in the foreseeable future. 
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Appendix A: Local Dynamics in LTB Spacetime 



Imagine a Keplerian binary system consisting of masses mi and m2 with positions Xi and X2 within the Fermi 
co-ordinate system associated with the LTB metric (57). The "Newtonian" tidal potential 

Vn = ^ K,iT) {X^ + Y') + l K,{T) (Al) 

is quadratic in X, so that the corresponding tidal acceleration, — VVn, is linear in X. This implies that the external 
cosmological perturbation on the Kepler system has the form 

F = -K^X + {K^ - K^) Z Z , (A2) 

where X = Xi — X2 describes relative position. Here, Ki = Ki{Tq) and K2 = K2{Tq) are cvahiatcd at the present 
epoch To. It is important to emphasize that our task here is to study the observable consequences of cosmological 
perturbations on local systems; therefore, the external perturbation on the binary system is evaluated at the present 
epoch in equation (A2). The equation of relative motion now takes the standard form 

d'X GMqX _ 

-^ + —^-^^ (A3) 

where Mq = mi +m2. We assume that the binary system experiences small perturbations due to the tidal influence of 
the cosmological gravitational field. The perturbed Kepler system (A3) can thus be treated using standard methods 
of celestial mechanics [23] . In our approach, the study of the long-term cosmological evolution of the Keplerian binary 
should be based upon the solution of (A3) in the more general case in which Ki{T) and K2{T) in F are not restricted 
to the present epoch; however, such an analysis is beyond the scope of this paper. 

The state of relative motion in (A3) is given by the position and velocity at a given time T; on the other hand, 
one could employ instead the six orbital elements of the instantaneous osculating ellipse to specify the motion. 
The temporal evolution of these orbital elements are described by the Lagrange planetary equations [23], which are 
therefore equivalent to (A3). For the orbital parameters, it is useful to employ Delaunay's action-angle elements 



(Z, G, H, I, g, given by 



L = A^/'^, G = [GMoA{l-e^)f'^ , H = Gcosi, (A4) 
? = 1? — e sini? , g = argument of the pericentre , h = longitude of the ascending node. (A5) 

Here, A is the semimajor axis of the osculating ellipse, e is its eccentricity, i is the orbital inclination, ■(? is the eccentric 
anomaly, and I is the mean anomaly. Moreover, G is the magnitude of the orbital angular momentum vector and H 
is its Z-component. Along the osculating ellipse, the radial position p can be expressed in terms of the true anomaly 
ip and eccentric anomaly respectively, as follows 



A{1- e") 
1 + e cos ip 



A(l-ecosi?). (A6) 



In equation (A4) and throughout, we consider only positive square roots. To express the dynamical equations in 
terms of Delaunay's elements, it is necessary to decompose the perturbing acceleration F in terms of an orthonormal 
frame field adapted to the osculating ellipse. That is, we can write 

F = Fpp + FsS + Fnfi, (A7) 

in terms of its radial, sideways, and normal components. Here p = X / p is the radial unit vector, s = n x p, and n 
is the imit vector in the direction of the orbital angular momentum G; that is, G = Gn. Thus, p and s are in the 
instantaneous orbital plane of the osculating ellipse, while n is normal to it. 
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The equations of motion in terms of Delaunay's elements are then given by 
dL 

— = [Fp e sin (fi + {1 + e cos ifi)], (A8) 



dG 
dT 
dH 
df 

dT 



pFs, (A9) 
p [Fg cos i — Fn sin i cos {cp + g)] , (AlO) 
= wk H [^P (~^^ "I" <P + e cos^ (^) — Fg (2 + e cos <p) sin (^] , (All) 



do pFn cos i . , _^ (1 — e^) / „ „ 2 + ecos(/? . , 

^ - ^" sin(vp + g)+ ^ / ( -Fp cosyj + fg ^ , __^ siny>), (A12) 

(A13) 



dT G sini ujYiAe \ l + ecos<^ 

dh _ pFn sin {ip + g) 
dT ~ G sini 



where wk = [GMq)^^'^ jl?' . A direct derivation of these equations is essentially contained in Appendix B of [24]. 
To find the explicit form of Fp, Fg, and F„ in terms of Delaunay's elements, we note that [24] 

X = p cos h cos {ip + g) — sin h cos i sin (ifi + g) , (A-14) 

Y = p sin /i cos {'fi + g)+ cos h cos i sin (f + g) , (A-15) 

Z = p sin i sin + g) . (A-16) 

Moreover [24], 

h = ^sin/i sini, — cos/i sinijcosi^ , (A17) 
s = ^— cos h sin {^p + g) — sin h cos i cos (v + ff) , 

— sin h sin ((^ + g) + cos /i cos i cos (<y5 + g) , sin i cos + g)^ . (A-18) 

A straightforward calculation using (A2) reveals that 

Fp = -K° p + {K° - K°) p sin^ i sin^ {(p + g) , (A19) 
Fs = {K^ - K^) p sin^ i sin ((^ + g) cos + g) , (A20) 
F„ = (fsTj' - if^) psini cosi sin((/7 + g) . (A21) 

Substituting these results in equations (A8)-(A13) and averaging over the "fast" orbital motion, one can determine 
the "slow" evolution of the orbit under cosmological perturbations in this case. 

The weak external perturbation (A2) naturally splits into two tcirms: one that is proportional to Ki and the other 
proportional to — The influence of these on the orbit can be analysed separately and the results can be 
superimposed in accordance with our linear perturbation scheme. Let us therefore first consider a tidal perturbation 
of the form —K^X. It turns out that this problem has already been solved in a different context [25]; in fact, 
the details of the averaging procedure are given in section IV of [25] for a perturbing acceleration of the form \X, 
where A is a constant. It is interesting to note that in [25], A = Ac^/3, where A is the cosmological constant 
associated with the Kerr-de Sitter spacetime. It is shown in [25] that the orbit is planar and its semimajor axis and 

1 /2 

eccentricity do not change on the average; however, the pericentre precesses with frequency 3A (l — e^) n/ (2wk). 
Thus the net average effect of the first term in (A2) on the orbit is to generate a pericentre precession of frequency 

-3/f?(l-e2)'/'n/(2u;K). 

The second term in (A2) is a perturbing acceleration in the Z-direction, which corresponds to the radial direction 
in the standard form of the LTB spacetime. Inspection of equations (A8)-(A13) reveals that dH /dT = 0, hence the 
Z-componcnt of the orbital angular momentum remains unchanged; moreover, there is clearly no effect to linear order 
if the orbit lies in the {X, Y) plane. Hence we assume i 0. It is then straightforward to average the right-hand 
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sides of the remaining equations over the "fast" motion with frequency ojk = 27r/TK such that for a function J^, 
{T) = Tj^^ T dt, or equivalently, 

^0 (l + ecos<p) 

using the unperturbed orbit. The resulting integrals can be evaluated in principle, but it is simpler to express the 
main results for a slightly eccentric orbit. We find that 

(^:i)=0{e^), (^)=Oie), (A23) 



,dT/ ^ ^ ' \dr, 

Jo \ (k" - K") 

^ ^ ^ ^ [cos^ i + (cos 2g-3 sin^ g) sin^ i] + O (e) , (A24) 



^ dT / 2ujk 
' dJi\ _ {K^-K°) cos^ 
dT / ~ 2wk 



+ O (e") . (A25) 



Thus the main effects here are the precessions of the pericentre and the ascending node with frequencies given by 
(A24)-(A25). The net average effect of (A2) on the orbit is mainly the precession of the orbit in its own plane as well 
as about the Z-axis (i.e. the LTB radial direction); the latter motion occurs with frequency (^K^ — K2) cosi/ (2cjk)- 
Finally, let us restrict our general off-centre LTB treatment to a situation where the fundamental observer is 
very close to the centre of spherical symmetry. Indeed, the maximum anisotropy in the temperature of the cosmic 
background radiation is at 0.002 level due to our motion relative to this radiation bath; therefore, within the LTB model 
a reasonable estimate may be obtained by assuming r = as in section V. Indeed, equation (62) implies that K1—K2 
for r — > 0. It follows that the cosmological perturbations on average leave the shape of the orbit and the orientation 

1 /2 

of the orbital plane unchanged, but cause a pericentre precession of frequency — 3go Hq (l — e^) n/ {2uj]^). For the 
motion of the Earth around the Smi, Hq/luk ~ 10^"^^ /yr, which is some fourteen orders of magnitude smaller than 
the Einstein pericentre precession. Thus the cosmological pericentre precession does not appear to be detectable in 
the foreseeable future. This conclusion applies to all of the cosmological models considered in this paper, so long as 
any deviation from isotropy about the Keplerian system can be neglected. 
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